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A method is presented for computing the compressive buckling load of a simply supported 
elastic rectangular plate having a central circular hole reinforced by a circular doubler 
plate. Numerical results are presented for six square plates having hole diameters up to 
one-half of the plate length. Comparison of these results with those computed for plates 
without holes shows that an unreinforced circular hole causes a relatively small reduction in 
buckling load, and reinforcement of a circular hole by a doubler plate causes a substantial 
increase in buckling load. 



I. Introduction 

The stressed skin cover of airplane wings and 
fuselages has to be perforated by holes at a number 
of places to give access to the interior of the wing 
or fuselage. These holes must be reinforced to 
prevent weakening of the entire structure by the 
hole. 

An ideal reinforcement would be one that 
confines the disturbance in the stress flow to the 
immediate neighborhood of the hole and which, at 
the same time, adds minimum weight to the 
structure. 

Too little is known about reinforcements to 
approximate this ideal in practice. In the 
absence of a rational design procedure, it has been 
customary to reinforce holes by circular doubler 
plates riveted to one side of the sheet, thereby 
increasing the effective thickness and rigidity of 
the sheet at the ex\^ of the hole. 

The Bureau of Aeronautics, Navy Department, 
initiated a study of reinforcements around holes 
at the National Bureau of Standards to provide a 
better understanding of reinforcements and to 
indicate practical improvements in their design. 
The first phase of this investigation was a plane 
stress analysis of a doubler plate reinforcement 



around a circular hole in an infinite sheet under 
uniform compression in one direction or under 
shear. 

The second phase of this investigation consisted 
of a check of the theoretical analysis by tests on 
plates having reinforced circular holes. It was 
found that the analysis gave an adequate descrip- 
tion of the stresses and displacements provided 
that the doubler plate was fastened to the sheet 
by at least two concentric rows of rivets. 

The third phase of this investigation is described 
in this paper. An analysis is given for the stabiliz- 
ing effect of the doubler plate for square plates 
with central circular holes when the plates are 
subjected to edge compression in one direction. 

II. Method of Analysis 

An energy method for determining the buckling 
load of rectangular plates of constant thickness 
under compressive loads is presented by Timo- 
shenko. 1 A review of this derivation shows that 
the method is also applicable to plates of variable 
thickness. In this case, Tinioshenko's integrals 
Ii and I 2 are 

i S. Timoshenko, Theory of elastic stability pp. 307 to 326 (McGraw-Hill 
Book Co., New York, N. Y., 1936). 
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where 

A = plate thickness (function of x and y). 
x, y= rectangular coordinates with origin at 
center of plate and x-axis in direction 
of load. 
w= lateral deflection of plate. 
Z}=M 3 /12(1- M 2 ), flexural rigidity of plate 

(function of x and y). 
jii = 0.3, Poisson's ratio. 
a x = tensile stress in x direction. 
(j y = tensile stress in y direction. 
t xv = shear stress. 
S= tensile stress in x direction far from hole. 

The stress ratios <r x /S, <r y /S, and r xy jS just prior 
to buckling may be obtained as shown by Gur- 
ney. 2 The stresses a r , a d , and r rd referred to 
polar coordinates r, are given, with some change 
in notation, as 



a r =:S{F+KR 2 /r 2 - (A+3CR i /r i +2DR 2 /r 2 ) cos 26' 
a 6 =S{F-KR 2 /r 2 + (A+6Br 2 /R 2 + ZCR*/r±) cos20] 
T r9 =S{(A+Wr 2 /R 2 -3CR 4 /r*-DR 2 /r 2 ) sin 20} 



where 



(2) 



r, 0= polar coordinates with origin at center of 
hole, 
0— angle between radius r to point and direc- 
tion of load, 
R= radius of hole. 

F, K, A, B, C, D= coefficients with different 
values in sheet and in 
doubler plate. 
The values of F, K, A, B, C, and D may be deter- 
mined from the dimensions for a particular 
doubler plate by using figure 1. The coefficients 
from "1" region, figure 1, are used in eq 2 when 
computing the stresses in the unreinforced sheet; 
the coefficients from "2" region are correspond- 

2 G. Gurney, Brit. Rep. Memo. No. 1834 (Feb. 1938). 



ingly used when computing stresses in the rein- 
forced sheet. The abscissa in figure 1, (T/t-1) 
(b 2 /a 2 -l), is the ratio of the volume of material in 
the reinforcement to the volume of material 
removed from the sheet to make the hole. 

The stresses in rectangular coordinates may be 
computed from the stresses in polar coordinates 
as given in eq 2 by using the conversion formulas : 



-a T cos 2 #+o0 sin 2 0— T r6 sin 



20 
<r v —<r r sin 2 0+o-0 cos 2 0+r r sin 20 

T zy =( * f sin 20 + Tro cos 26. 



(3) 



The lateral deflection w of a simply supported 
rectangular plate of length a and width b will be 
approximated by the first terms in the trigono- 
metric series 



TX TV . TX ?>Tlj . 

w=a n cos — cos -r- + a 13 cos — cos — jf-+ 
a o a o 

3wx tv . 3tx Ztv . 

a n cos — cos -7-+a 33 cos cos -i-^+ . 

a b a b 



(4) 



Substitution of eq 4 in eq 1 and integration 
(see next section for method of numerical in- 
tegration) reduces integrals I x and 1 2 to quadratic 

expressions in the coefficients a n , #13 

The critical value of S at which buckling of the 
plate occurs is obtained from these expressions, 
according to Timoshenko, as that value of S 
which reduces to zero the determinant of the 
coefficients of a n , a^, etc. in the set of simultaneous 
equations, 



oa n oa n 
ctoi3 d«i3 



(5) 



The sign before S in eq 5 is plus, and that before y 
in the corresponding eq 212 given by Timoshenko 3 
is minus. Where S is considered positive for 
tensile stress, as is done in this report, the signs in 
eq 5 are correct. 

3 S. Timoshenko, Theory of elastic stability (McGraw-Hill Book Co., 
New York, N. Y., 1936). 
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Figure 1. — Values of coefficients F, K, A, B, C, and 1) to be used in eq 2 for computing the stress distribution in the 

sheet and reinforcement. 
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III. Numerical Integration 

The evaluation of the integrals in eq 1 over the 
surface of the plate is made difficult by having a 
circular inner and a rectangular outer boundary 
and by involving a stress that is a complicated 
function when expressed in rectangular coordi- 
nates. 

The integral 2Y, eq 1, was evaluated as the 
resultant of three integrations : 



ll — J-\a\J-lb Ac 



where 



(6) 



7 la = The integral I x for the surface enclosed 
by the outer rectangular boundary of the 
plate taking D=D iy the value of D in 
the unreinforced portion of the plate. 

jTi & =The integral I x for the surface enclosed 
by the outer circular boundary of the 
reinforcement taking D=D 2 —D U where 
D 2 is the value of D in the reinforcement. 

/i c =Tlie integral I Y for the surface enclosed 
by the circular boundary of the hole 
taking D=D 2 . 

The double integration necessary to evaluate I la 
could in every case be done directly. However, 
this was not possible in evaluating I lh and I u . 
The integrals I lb and I u were obtained by integrat- 
ing directly in respect to one variable and using 
Gauss' method of numerical integration, 4 for in- 
tegration with respect to the other variable. In 
each case where numerical integration was used to 
evaluate T u the number of Gauss points was in- 
creased to the place where the addition of another 
point caused less than 1 percent change in the 
integral. In general, five Gauss points were 
enough for this purpose. 

The integral I 2 was evaluated by considering 
the surface of the plate in three portions. The 
first portion, A, figure 2, was taken as the circular 
disk of the reinforced area; the second portion, 

B, figure 2, was taken as the circular disk between 
the outer edge of the reinforcement and the largest 
inscribed circle in the plate; and the third portion, 

C, figure 2, was taken as the remainder of the 
plate. The integral for each circular portion was 
determined, using Gauss' method of numerical 
integration, by first integrating numerically in a 

4 E. T. Whittaker and G. Robinson, Calculus of observations, p. 159 
(Blackie and Son, Ltd., London, 1944). 




Figure 2. — Subdivision of plate into three portions A, B, C, 
in evaluating I 2 . 

circumferential direction and then integrating 
numerically in a radial direction. The integral 
for the remainer of the plate was obtained also by 
using Gauss' method, first integrating numerically 
in the direction of the load and then at right angles 
to the direction of the load. Only one quadrant 
of the plate had to be considered because of sym- 
metry. In each case, a sufficient number of 
Gauss points was used to reduce the estimated 
error to less than 1 percent. Twenty- three points 
in*each quadrant were used for cases 1, 3, 4, and 
5, table 1, thirty-three points for case 2, and 
twenty-eight points for case 6. 

Table 1. — Dimensions for plates investigated 



Case 


Width 


Radius of 
hole, R 


Outer 
radius of 
reinforce- 
ment 


Thickness 


Direc- 
tion of 
load 


Trans- 
verse to 
load 


Rein- 
forced 
region 


Unrein- 
forced 
region 


1 


a 


a 


0. 125a 


0) 


0) 


h 


2 


a 


a 


. 0625a 


0) 


(i) 


h 


3 


a 


a 


.25a 


0) 


(») 


h 


4 


a 


a 


.125a 


0.25a 


2. Oh 


h 


5 


a 


a 


.125a 


.25a 


Lfift 


h 


6 


a 


a 


.125a 


. 1875a 


3.4ft 


h 



1 Hole in plate not reinforced. 

A typical distribution of Gauss points for 
computing I 2 is shown in figure 3 for one quadrant 
of a square plate with a central hole. For this 
distribution, the integral of a function F over the 
plate surface is 
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Figure 3. — Distribution of Gauss points for computing 7 2 in case 4- 
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/ - /'jp , dA=a 2 [0.075316jPi + 0.068366F 2 + 
surface 

0.068366^3+0.041605F 4 + 

0.084048F 5 +0.099898^o+ 

0.084048F 7 +0.041605^ 8 + 

0.028176^ 9 +0.056919^o+ 

0.067653^ + 0.056919^2+ 

0.028176F 13 + 0.010401^ 14 + 

0.021012^5+0.024975^16+ 

0.021012i< 7 17 +0.010401i<\ s + 

0.007044jF 19 +0.014230^ 20 + 

0.016913^2i + 0.014230F 22 + 

0.007044^ 23 ] 



(7) 



where Fi, F 2 . . . F 2 % are the values of the 
function F at the positions indicated in figure 3. 
An approximate check on the over-all adequacy 
of the numerical integration methods used was 
obtained by determining the buckling load of a 
square plate with no hole both by the exact 
method, 5 and by numerical integration with 23 
points for I 2 . The results differed by only 
0.7 percent. In the case of a plate with a rein- 
forced hole, it is likley that the more complicated 
stress distribution and the greater prominence 
of higher order terms in the series used for the 
deflection causes the error to be somewhat higher. 

IV. Convergence of Deflection Function 

The correct value of S for buckling of the plate 
would be that value of S which reduces the deter- 
minant of the coefficients of a n , a 1Z:> ... in the 
infinite set of eq 5 to zero. In order to limit the 
work of computing Sto a finite amount, preliminary 
computations were made to see which terms in 
the deflection function, eq 4, were most important 
in determining the value of S for buckling. These 
computations were made for the square plate of 
case 4, table 1. The compressive stress (-S) was 
10.33 Eh 2 /a 2 using only the first term of eq 4, 



w= 



TX TV 

--a n cos - — cos -T* 
a o 



Using only terms 1 and 2 of eq 4, 



TX TV 

w=a n cos — cos -^ 
a 



TTX 3 TTV 

-a 13 cos — cos - 



b a b 

(-S) was 9.43 Eh 2 /a 2 , a decrease of 8.6 percent. 
Using only terms 1 and 3, 

6 S. Timoshenko, Theory of elastic stability, p. 331 (McGraw-Hill Book Co., 
New York, N. Y., 1936). 



TX TV . 3TX 7TV 

w=a n cos — cos -f- +d 31 cos cos -A 

(10 dO 

(-S) was 8.80 Eh 2 /a 2 , a decrease of 14.7 percent. 
Using only terms 1 and 4, 

TX TV . 3tx 3tv 

w= a X i cos — cos-r+a 33 cos cos^-^> 

a b a b 

(-S) was 10.32 Eh 2 1 a 2 , a decrease of 0.1 percent. 
Had additional terms of the trigonometric series 
been used with term 1, the decrease in the buckling 
stress would probably have been proportionately 
smaller. It is believed, on this basis, that the 
first three terms of eq 4 approximated the lateral 
deflection of the square plates with simply sup- 
ported edges with sufficient accuracy to give the 
buckling stresses within 5 percent. Accordingly, 
only the terms 

tx Ty . TX 3wy . 

W = a n COS COS ^ + #13 cos — cos — r^ + 

a b a b 

3tx Ty 
a 3 i cos cos --r- 

a b 

were used in the remaining computations for 
square plates. 

V. Results 

The compressive buckling load was computed 
for six square plates with reinforced and unre- 
inforced holes, figure 4. The dimensions for the 
plates are given in table 1. 

The analysis gave the compressive buckling 
stress (—8) far from the hole corresponding to 
the stress distribution derived in reference [2]. 
Values of (— S) are given in table 2. The com- 
pressive stress (-S) is somewhat larger than the 
average stress in the direction of the load along 
the edge of the plate. A better value of the 
average buckling stress <r cr along the edge of the 
plate was obtained from 



(71 + 0-2 



(8) 



where <j x is the average compressive stress on the 
loaded edge of the plate and <r 2 the average com- 
pressive stress on the center cross section, obtained 
by dividing the load on that section by ah. The 
value of <7i was obtained by numerical integration 
and the value of a 2 by direct integration of the 
stresses a x given by eq 3. Values of (-S), a u <r 2 , 
and <r cr are given in table 2. The average stresses 
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Ficure 4. — Loading of square plates with simply supported 
edges. 



Dimensions in table 1. 



o*! and <r 2 differed from each other by less than 3 
percent, indicating that the equilibrium condition 
is substantially satisfied. 

Table 2. — Critical compressive stress for plates investigated 



Case 


Stress (-.S) 
far from hole 


A. verage stress 
<r\ at 

. a a a 
/ %'~2° < 2 


Average stress 
o"2 correspond- 
ing to load at 
i=() divided 
by ah 


Buckling stress, 

<r C r, average of 

<n and a. 


1 
2 
3 
4 
5 
6 


3.611 £/j2/o 2 
3.620 Ehya* 
3.750 EM/at 
8.021 £/i 2 /« 2 
5.700 Ehya* 
7.945 EhVa* 


3.454 E/t 2 /a2 
3.578 Eh^/a* 
3.076 Eh 2 /a 2 
7.984 M 2 /a2 
5.578 M 2 /a 2 
7.881 Eh^a 2 


3.490 EftVo 1 
3.591 EM/a* 
3.164 ^2^2 
8.069 E/i 2 /a2 
5.646 MVa 2 
7.984 M2/a2 


3.472 £/i2/a 2 
3.584 E/(2/o 2 
3.120 Ma/oa 
8.027 M»/b> 
5.612 E/i2/a 2 
7.932 M»/a 2 



The buckling stress of a simply supported 
square plate of constant thickness without a hole 
is given on p. 331 of reference [3] (taking Poisson's 
ratio as 0.3) as 

<r cr =3.615 Eh 2 ja 2 . (9) 



Comparison of this value with the values computed 
for cases 1 to 3 of table 2 shows that unreinforced 
central holes, having a diameter of Y%, % and )i of 
the length of a side of a square plate, reduce its 
buckling load by only 1, 4, and 14 percent, 
respectively. Comparison of the value in eq 9 
with the values for cases 4 to 6 of table 2 shows 
that the reinforcements increased the buckling 
stress over that for the plate without a hole by 55 
to 122 percent. 

The effect of the shape of the reinforcement on 
the buckling load is indicated by comparing cases 
4 and 6, table 2. For these cases, the hole size and 
volume of reinforcing material are identical but 
the reinforcing material is concentrated nearer the 
edge of the hole in case 6 (see fig. 4). The com- 
puted buckling stresses are the same within 1.2 
percent. 

VI. Conclusions 

A numerical procedure was developed for 
estimating the buckling stress of simply supported 
rectangular plates with circular holes and doubler 
plate reinforcement. The procedure provides a 
convenient method for solving the integrals for 
the energy stored in the plate. 

The computations showed that the buckling 
stress of square plates is reduced only a small 
amount by the presence of unreinforced holes. For 
the cases considered, the greatest reduction, 14 
percent, corresponded to case 3, a square plate with 
a hole diameter 0.5 times the length of the plate. 

Reinforcements of the doubler-plate type are 
shown by the computations to cause marked 
increases in the buckling stress. For case 4 with 
reinforcement, the computed buckling stress was 
2.3 times that for case 1 with no reinforcement. 

The buckling stress seems to be insensitive to 
change in shape of reinforcements of given 
volume. The computed buckling stresses were 
the same within 1.2 percent for cases 4 and 6, 
which differed in having doubler plates of different 
thickness and radius, but of the same volume. 



The authors thank the Bureau of Aeronautics, 
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to publish the results. 
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